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RESUMEN 

Aquf presentamos soluciones semianaliticas de la hidrodinamica de remanentes de supernova compactos 
(cSNRs) y su evolucion qufmica. Discutimos dos modelos de explosiones de supernova de 15 M y 25 Mq, 



envueltas en un medio denso (no 



10' 



'). Demostramos que despues de algunos anos de evolucion, la 



difusion produce una distribution homogenea de los productos de la SN dentro de los delgados cascaroncs 
frios formados despues del enfriamiento del material barrido. Encontramos metalicidades en un intervalo de 
5.5 < Zo/Zq < 6.0 para M$n = 25 Mg y 1.0 < Zo/Zq < 1.5 para Msn = 15 Mg, y especulamos que 
cstas metalicidades estan relacionadas con las abundancias observadas en QSOs (Fcrland et al. 1996, Hamann 
& Fcrland 1999). 

ABSTRACT 

Here we present semi-analytical solutions of the hydrodynamics of compact supernova remnants (cSNRs) and 
their chemical evolution. Two models of 15 Mq and 25 Mq supernova explosions evolving in a high density 



medium (no ~ 10' 



are throughly discussed. We demostrate, that after a few years of evolution, diffusion 



leads to an homogeneous distribution of the SN products inside the cold, thin shell formed after cooling of the 



swept-up matter sets in. Wc find mctallicity values within a range 5.5 < Zo/Zq < 6.0 for M, 
and 1.0 < Zq/Zq < 1.5 for M ; 



25 M, 



SN 



15 M , and speculate wheather these mctallicitics could be related to 



the observed QSOs super solar elemental abundances (Fcrland et al. 1996, Hamann & Fcrland 1999). 
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1. INTRODUCTION 

The discovery of the high-rcdshift quasars (QSO) 
implies that massive star formation occurs very early 
in galaxy evolution. Spectra of high rcdshiftcd QSOs 
(z > 4) present broad line emission (Hamann & Frc- 
land, 1993) which is similar to that seeing in active 
galactic nuclei (AGNs). The metallicitics of the as- 
sociated fast moving gas derived from the spectral 
analysis are surprisingly high for so young objects 
and fall into the range Z Q < Z < 9Z Q (Hamann &; 
Ferland 1999). These results have a great interest 
because provide information about an early history 
of star formation in galaxies. If high metallicity is 
related to a short episode of star formation, it has to 
occur at early (few Gyr) stages of galaxies evolution 
and suggests a fast mixing of ejected metals with a 
bulk of primordial low metalicity gas (Hamman & 
Ferland, 1993; Ferland et al. 1996). Different modi- 
fications of this scenario have been proposed by Silk 
& Rees (1998) who have discussed an inhenced star 
formation inside an accretion disk. 

The observed quasar (and AGN) spectra show 
a great similarity to the spectra of type II peculiar 
supernova remnants (Terlevich et al. 1992) which 
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one believes to be a compact SNR (cSNR) evolving 
in a very dense (~ 10 7 cm -3 ) circumstellar envi- 
ronment. Therefore here we analyze an alternative 
scenario based on the starburst model of AGNs (Ter- 
levich et al. 1992). For the calculations we use a 
semi-analitycal model of the young supernova rem- 
nant evolution and thoroughly discuss the dispersal 
of the ejected metals into the swept-up primordial 
gas. We simulate cSNRs evolution from the 15 Mq 
and 25 Mq progenitors and show that model pre- 
dicted thicknesses, densities and temperatures allow 
a fast mixing of the abundant ejecta with the swept- 
up low mctallicity circumstellar material. The resul- 
tant cSNR abundances fall into the observed QSOs 
metallicity range. 

2. THE INPUT MODEL 

2.1. Ejecta model 

The interaction of the SN ejected matter with the 
surrounding medium causes a four zone structure: 
the inner free expanding ejecta which is bounded by 
a reverse shock with radius Rrs and velocity D#s; 
the shocked ejecta which is accommodated between 
the reverse shock position and the contact discon- 
tinuity Ren; the shocked CSM which is separated 
from the ejected material by the CD and from the 
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surrounding medium by the leading shock with ra- 
dius Rls and velocity Dls- The outward shock 
sweeps-up the dense surrounding medium and en- 
hances its temperature to ~ 10 8 K. The reverse shock 
is slower and thermalizes the SN ejecta. 

The cSNR initial evolution is dependent on the 
surrounding gas number density and the ejecta den- 
sity and velocity radial distribution. We adopt 
Franco et al. (1991) analytical model based on 
Wooslcy ct al. (1988) calculations. It assumes that 
supernova ejecta obey the free expansion law with 
r~ 3 density and linear velocity distributions: 



Pej 0, t) 



V ej (r,t) 



r <R r 



(1) 



r-R c T/ 
R{t)-R a Vm > 





r>R c 
r < Rc 



(2) 

where R(t) = R m + V m t, R m is the initial radius 
and V m is the maximum ejecta velocity, R c is the 
inner radius of the ejecta. The total mass of ejecta 
and masses of the ejected oxygen and iron have been 
taken form Woosley (1978) model. The initial model 
parameters are summarized in the Table 1. 

TABLE 1 
INITIAL CONDITIONS 





15 M Q 


25 Mq 


Initial energy [erg] 


1 x 10 51 


1 x 10 51 


Ejecta mass [Mq] 


13.44 


23.61 


Maximum 


6000 


5000 


velocity [km s _1 ] 






R c [cm] 


6.45 x 10 15 


3.24 x 10 15 


Rm [cm] 


3.82 x 10 16 


4.10 x 10 16 


Ejecta mass 


0.42 


3.0 


of oxygen [Mq] 






Ejecta mass 


0.06 


0.08 


of iron [Mq] 







2.2. The dispersal of metals 

There are essentially two irreversible mass trans- 
port mechanisms leading to the dispersal of the 
heavy elements released during the supernova explo- 
sions throughout the ISM. They are diffusion and 
turbulent mixing (Roy & Kunth, 1995; Tcnorio- 
Tagle, 1996, Ocy 2003). The determination of the 



appropriate turbulent parameters so far is specula- 
tive. Therefore here we focus on the diffusion of the 
supernova ejected metals and their mixing with the 
swept-up CSM. 

We assume that the ejecta represents a homo- 
geneous mixture of the progenitor outer hydrogen 
envelope with the synthesized heavy elements. At 
the contact discontinuity this high metallicity gas is 
mixed with the low metallicity circumstellar mate- 
rial. The characteristic diffusion time over a length 
L is to = L 2 /D, where D is the diffusion coefficient. 
When the dynamical time t exceeds the characteris- 
tic diffusion time tu, the metallicity of the mixture 
approach the value 



Z A {t)/Z Q 



M A eJ (t)/Z A q + Zism M ISM {t) 
M ISM (t)+ X M ej (t) 



(3) 



where McsM(t) and %M ej -(i) are masses of the 
swept-up CSM and the ejecta, respectively, and 
M A ej (t) is the mass of the element A that has been 
accumulated within a given gaseous segment. 

3. THE RESULTS OF THE CALCULATIONS 
3.1. Hydrodynamical evolution of the cSNR 

Figure 1 shows runs of leading and reverse shocks 
velocities and radii for 15 M Q supernova remnant. 
The leading shock has initial speed ~ 4200 km s _1 . 
Strong radiative cooling in the swept-up CSM mat- 
ter induces a sudden drop in the post-shock thermal 
pressure and a noticeable deceleration of the lead- 
ing shock velocity which falls down from ~ 3000 
km s _1 to ~ 2000 km s _1 during a short transi- 
tion phase. The process is followed by the condensa- 
tion of the swept-up interstellar gas into a thin dense 
shell. The shell density exceeds the shocked ejecta 
density by about four orders of magnitude. There- 
fore the CD adjusts the leading shock position and 
later on follows its expansion. Radiative phase starts 
at t sg ~ 1.2yr after explosion and is well identified 
as breaks of the solid lines at velocity curve. The 
reverse shock is substantially slower. It remains adi- 
abatic for a much longer time because of the smaller 
post-shock density and thus less effective cooling. 
Shock waves heat both the swept-up circumstellar 
material and the shocked ejecta to high temperatures 
~ 10 8 K. However the temperature of the shocked in- 
terstellar material suddenly drops when cooling sets 
in and later on is supported at 10 4 K level by the ion- 
izing photons arising from the cooling shocked gas. 

Figure 2 shows the swept-up mass time evolution. 
The interstellar mass accumulated in the outer shell 
grows almost linearly with time and at the end of the 
calculations reaches around 10 Mq. During this time 
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Fig. 1. Radius and velocity evolution of the shocks. 
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Fig. 3. The cSNR metallicity as function of time if oxy- 
gen is used as a tracer. 



the reverse shock processes only ~ 4 M Q of the ejecta 
which is spread inside a relatively thick inner shell. 
The flat parts of the curves in Figure 2 come from 
our simplified description of the swept-up interstel- 
lar mass condensation when catastrophic radiative 
cooling sets in. 
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Fig. 2. Swept-up mass evolution. 

3.2. Chemical evolution of the cSNR 

Using the shell thickness and temperature given 
in the previous section we calculate the character- 
istic diffusion time t<j one needs to mix the ejected 
metals with the swept-up CSM. At the beginning 
of the calculations the characteristic diffusion time 
grows with time due to the growth of the shell thick- 
ness. However when radiative cooling sets in it drops 
about four orders of magnitude because of the sud- 
den decrease of the outer shell thickness. This leads 
to an effective and rapid mixing of the ejected met- 
als with the swept-up interstellar matter and defines 
the fast chemical evolution of the remnant. Initial 
mctallicitics of the swept-up interstellar medium and 
the ejecta and the resultant metallicity of the cSNR 
outer shell arc shown in Figure 3. 



4. CONCLUSIONS 

The calculations show that effective metal 
diffusion occurs at early cSNR evolution, soon 
after a dense shell formation. The resultant shell 
metallicity was found to be within the range 
1.0 < Z/Zq < 6.0 if oxygen is used as a tracer 
and 0.65 < Z/Zq < 1.3 if iron is used as a tracer 
of the gas metallicity. This is in a remarkably good 
coincidence with the observed QSO metal abun- 
dances. We speculate therefore that the observed 
QSO mctallicitics could be related to SN explosions 
immersed in a high density CSM. 
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